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We consider resonant transmission through a finite-length quantum wire connected to leads via 
finite transparency junctions. The coherent electron transport is strongly modified by the Coulomb 
interaction. The low-temperature current- voltage {IV) curves show step-like dependence on the bias 
voltage determined by the distance between the quantum levels inside the conductor, the pattern 
being dependent on the ratio between the charging energy and level spacing. If the system is 
tuned close to the resonance condition by the gate voltage, the low-voltage IV curve is Ohmic. At 
large Coulomb energy and low temperatures, the conductance is temperature-independent for any 
relationship between temperature, level spacing, and coupling between the wire and the leads. 
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Quantization of conductance in ballistic one dimen- 
sional (ID) channels or islands [1[ is one of the central 
issues of current mesoscopic physics, see ^ for a review. 
In small-size islands, the Coulomb blockade suppresses 
conductivity at low temperatures and small applied volt- 
ages |3|. The Coulomb blockade can be controlled by 
varying the charge of the island between two tunnel bar- 
riers by the potential at the gate electrode; the conduc- 
tance of such system shows oscillations as a function of 
the gate voltage, due to the periodic modulation of the 
charging energy. The classical theory of the Coulomb- 
blockade oscillations was developed in [j|, while the role 
of the discreteness of the spectrum was addressed in [5|, |6[ 
and several subsequent papers. In these models the is- 
land was assumed to be almost isolated from both source 
and drain, so that the number of particles on it could 
be considered a conserved quantity. This assumption, in 
general, does not hold if the transparency of the contacts 
between the island and the leads is finite Ml. Typical 
examples of such systems are quantum dots [8| . 

Transport through quantum conductors (QCs) essen- 
tially coupled to the leads has not been fully investi- 
gated so far, though it offers excellent opportunities for 
studying the interplay between quantum and classical 
properties of QC. In this Letter, we investigate trans- 
port through a relatively short single-mode QC assuming 
a simplest model where charging effects are important 
whereas the Luttinger-liquid behavior is still not pro- 
nounced (see, e.g., [9[). We consider a ballistic conductor 
connected to two leads via identical contacts with a fi- 
nite transition amplitude T. Its length d is shorter that 
the electron mean free path, the transport mechanism be- 
ing the resonant transmission, which we distinguish from 
elastic cotunneling [lO|. This mechanism is relevant to 
recent experiments on carbon nanotubes [ll|, \l^ . 

At T -^ 0, the QC holds a fixed number of particles. 



N. If a particle with energy E tunnels between a lead 
and the state with energy ep in the QC, the charging 
energy of the QC couples E and ep through [6[ 

E^ep + 2Ec{N-No)±Ec (1) 

where ± stands for adding/removing a particle, Ec = 
e^/2C is the Coulomb energy of the QC, C is the QC 
capacitance, and eA^o is the gate-induced charge density. 
At finite T, the system forms a double-barrier resonant 
tunneling structure, its properties being determined by 
the resonant states. The distance between these resonant 
levels is AE = Trhv/d (where v is the electron velocity) 
while their width is F = \T\'^AE. The occupation prob- 
ability of the modes is determined by the competition 
between the relaxation processes in the QC and in the 
leads, as well as by escape probability from the QC into 
the leads. If the inverse escape time, T/h, is much less 
than the inter-mode relaxation time, t~^ , the occupation 
probability is the equilibrium Fermi distribution with the 
lattice temperature, the chemical potential being deter- 
mined by the (conserved) number of particles. In the 
opposite case where F ^ h/r the distribution inside the 
QC is determined by the coupling to the leads; it can be 
strongly nonequilibrium. The distribution function is to 
be found from the kinetic equation. The results, in gen- 
eral, depend on interplay between the different relaxation 
mechanisms in the quantum conductor and the leads. 

We will assume that the resonance width F is small, 
but is still larger than the relaxation rate in the QC, 



?i/T «: F < A£; . 



(2) 



The requirement ^ is opposite to the condition F ^ h/r 
considered in [6|. In this limit QC cannot be treated as 
an isolated quantum dot with the vanishingly weak cou- 
pling to leads [lj|, and therefore, A^ is not a good quan- 
tum number but rather an average value determined by 



the interaction with the leads. To reflect this we employ 
the Hartree-Fock type approach which leads to modifica- 
tion of Eq. ([T]). We show that at finite T the excitation 
spectrum changes significantly. At large Coulomb en- 
ergy, Ec > Ai?, and low temperatures, T <C Eq, energy 
exhibits a sharp step as a function of the internal mo- 
mentum in the QC. This step defines the internal Fermi 
level. The width of the step is determined by the width of 
the resonant level. As a result, zero-voltage conductance 
becomes temperature-independent already at T ^ Ec 
irrespective of the relationship between T, Ai?, and F. 
In a carbon nanotube placed between the two leads the 
ratio Eq/^E = e^d/2TThvpC can reach values e^ /2iThvp 
for the minimal capacitance of the tube C '^ d. For typ- 
ical vp = 0.8 • 10^ cm/s 1^ this ratio is ~ 0.46, i.e., is of 



the order of unity. Therefore, in order to achieve better 
understanding of the experimental situation and having 
in mind more general applications, both limits of large 
and small ratio Ec/ ^E have to be studied. 

Model.- We specify the charging Hamiltonian as 



He — Er 



J2 f i'i{rH^ir)dV-No 
„. Jv 



(3) 



V is volume of the QC and a is the spin index. The spin 
dependence is due to the level filling that controls the 
charge on the QC. The spectral and transport properties 
are determined by the retarded (advanced) and Keldysh 
Green functions satisfying the equations 
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= G^\e,r,) 
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Gf(ri,r2) 



Go '(e, ri) = -e + r|2p2/2m - ^ + f^C Pi = -id/dn. 

The energy Uc — e'^{N — Nq)/C is produced by the 
charge on the QC, /i is the chemical potential; ri belongs 
to the QC, otherwise the charging energy vanishes. 

For methodical purposes, let us start with the example 
of infinitely long QC assuming AE <^ Ec- Since Ec is 
coordinate independent one can consider the wave func- 
tions as plane waves and use momentum representation. 
Then the Green function has the pole at e = Cp where 

Sp = (p + Uc+Ectaiih^ , ^p = ^-Ef . (4) 

Ep is measured from the Fermi level /i = Ep in the leads. 
At Ep 3> T this coincides with Eq. H]). For Cp — > one 
has Sp = {^p + Uc)/{l-Ec/2T). The slope of the Sp vs Cp 
dependence is negative for Ec > 2T, and the function Sp 
becomes multi- valued, see Fig.[TJa). As we will see later, 
Ep in fact experiences an abrupt jump as a function of 
internal momentum shown in Fig. [ija) by the vertical 
line. This infinite slope is a consequence of the chosen 
approximation of an infinitely long QC with AE — > 0. 



In a finite QC, the jump acquires a finite width de- 
termined by the transparency of the contacts. To study 
this problem in more detail it is convenient to expand 
Green's functions over the orthonormal set of functions 
M„(r), which satisfy the Schrodinger equation 



2m 



l^ + Uc 



j Un + £^C ^ UmNmn = EnUn 



(5) 



Here A^„„ = /dV<„(r)«„(r)/(£;™); f{E^) = 1 - 
2n[Em), where n{Em) is the occupation probability of 
the m-th state. The diagonal element is proportional to 
the average charge eA^„„ in the state n. 

Spectrum. - Since main contributions to the trans- 
port come from the quasi-resonant tunneling states it is 
convenient to expand the functions u„ over the scattering 
plane- wave states, Wp(a;), satisfying the equation 



2^2 



2m 



V^5{x±d/2) 



^p{^) = ^pVpix) . (6) 



The wave functions can be chosen as incident waves on 
the left, Vp, and incident waves on the right, Vp. We 
assume a symmetric structure such that each barrier is 
characterized by the same plane-wave reflection TZ and 
transmission T amplitudes, TZ = \TZ\e^^ and T = — |T|e** 
where S is the scattering phase; \TZ\'^ + |Tp = 1. We now 
express the wave functions u„ through the free functions 
Vp: Un{x) = {2Tr)~^ J AnpVp{x) dp. Using the orthogonal- 
ity and completeness of both sets m„ and Vp one can show 
that the states u„ with the different n are expanded in 
Vp with different p. This property, in particular, excludes 
the multi- valued solutions of Eq. (j4|) shown in Fig.[lja). 
Equation ([5|) can be then rewritten in the form 

Cp - En= / T^-r^ [(e</? - Uc)Mpp, - EcQpp'] ; 

J 211 Anp 



rd/2 
Mpp' = / v*{x)vp'{x)dx, 

J-d/2 



(7) 



^PP' — / y / 9 O ^nipi^mp2J lv^mj-'"PPl-'"P2P'■ 
The matrix elements Afpp' can be explicitly expressed 
through the transition amplitude, T. 

In the equilibrium, the distribution functions for par- 
ticles coming from the left and right are the same, 
f{E) = ta,nh.{E/2T). In what follows we concentrate 
on the situation when the transparency satisfies Eq. ^. 
For \T\ <^ 1 there exist sharp resonances in the trans- 
mission when the momentum of electrons inside the QC 
is close to the resonant values corresponding to integer 
ratio s = {pd + 5)/tt, 5 k, it. We denote p^ the even and 
odd resonant states for s = 2?^ + 1 and s = 2n, respec- 
tively. At Ec ^ AE these internal resonance momenta 
of w„(r) lie far from the incident particle momenta in the 




FIG. 1: (Color online) (a) Energy spectrum of an infinitely 
long normal conductor. Finite temperature; The slope at 
_B = is positive for Ec < 2T (green line) and negative for 
2T < Ec (dotted line). For T = (red lines) the Fermi sur- 
face corresponds to the jump in the spectrum; pj^ = Pp. — Uc /v 
(dashed vertical line). Straight (black) line: zero Coulomb en- 
ergy, (b) Energy spectrum for a finite-transparency contacts. 
The dots illustrate sequential position of even and odd states. 



leads y/2m{E -\- Ep) but arc close to the resonance mo- 
menta of the plane- wave states, Vp{x). This implies that 
transport is controlled only by discrete resonant levels, 
which distinguishes the considered mechanism from the 
usual elastic cotunneling through a multi-level quantum 
dot [10| . Consequently, one can put A 



np 



M* 
gr^p 



KL 



pqn 



where the continuous parameter qn is the wave vector of a 
particle inside the QC close to n-th resonance. Equation 
([7|) for p = g„ becomes (for both even and odd states) 

dp 



■Uc - e(p + Ec 



M. 



■'p ■'p 



27r 

gi{p-'P')d/2^lp 



fiEp)M, 



± . 
pp ■ 



(8) 



pp 



(g—ip'd/2 q; 'J^gip'd/2\/^ipd/2 ^ J^* g — ipd/2\ 
-2^2 



Here €q^ = h qf^/2m — Ep, the integration is performed 
over the vicinity of the n-th resonance, p„; A(p) = 
(2/p) sin(pd/2). Equation ([8|) generalizes Eq. ((4|) with 
the replacements Sp — > Em, S,p — > £q„ and defines Qn cor- 
responding to the resonant transmission through the QC. 

If the intra-QC relaxation is efficient, t^ <^ ?i/r, where 
Tg is the inelastic relaxation time in the QC, the distri- 
bution is determined by the states in the QC, f{Ep) = 
fi^qn) = tanh(eq„/2T). The abrupt dependence shown 
in Fig. [DJa) is then smeared by the finite temperature, 
and En as a function of Qn acquires a finite width '^ T, 
similarly to the situation considered in [6|. However, if 
Eq. ^ holds, and when temperature satisfies T <^ Ec, 
the spectrum En{qn) shows a jump over 2Ec at some 
value qn = p^ such that all the levels with g„ > p^ are 
empty, while the levels with g„ < p^ are occupied. There- 
fore f{En) ~ sign((7„ —pp). With this approximation in 
the r.h.s. of Eq. ^ one finds for T < 1 

En = fiv{qn - pp) -eip + Uc + Ec^Pn - P^) ■ (9) 

Here g„ is close to the resonance value p„ = im/d and 

^{Pn ~ Pm) = (2/^) arctan [2d(p„ - P^)/|r|2] . (10) 



Equation ([9]) generalizes Eq. ((T|) for low transparency lim- 
ited by Eq. (I2|); Eqs. ([9]) and dJ) coincide far from the 
Fermi level. Energies En for the resonance momenta Pn 
are shown in Fig. [Ijb). The width of this function is 
determined by F, rather than by the temperature. 

The Fermi momentum, p^, is related to the number of 
electrons inside the QC. We define the variation of the 
particle number, SN , as a function of the bias and gate 
voltage. Using Eqs. ([9]), (|10p one can cast this variation 

in the form 5N = -{1/2)6 J2n a *(?'" ~ Pm) ■ Here the 
sum runs also over the electron spin index a. Let us de- 
note the average number of electrons for zero gate voltage 
as iVo. If A^o is an integer, the chemical potential lies far 
from one of the resonances \T\'^ /d ^ |p^ — p„| < Trn/2d. 
For non-integer A^Oi one can split it into the integer part, 
[iVo], and the remainder, iVo. In this case, the chemical 
potential is close to one of the resonances. Assuming that 
the state n is the closest to the Fermi level one can deter- 
mine p^ from the equation $(pn ~P^J.) — 1 ^ 2No. In gen- 
eral, the chemical potential and the number of particles 
are related by the equation 2{N — [No]) + ^{pn —Pfi) = 1. 
Near the degeneracy point when N — [Nq] = ±1/2, the 
Fermi momentum is close to one of the resonances. 

Current. For a finite bias V, energies in the left and 
right leads are shifted by ±6^/2. Following the same 
approximation as above, we have instead of Eq. Q 



En = hvp{qn - pp) - eip + Uc + 
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Now the Fermi momenta, p^ and p^, are not equal but 
are determined by the states in the left and right elec- 
trode, respectively. Similarly to the equilibrium case, we 
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Defining 



jp = u*{x)[dup{x) / dx] one can express current as 

^ = -^ E [A'^(^p) I'^-?P + ^''(^p) Imjp] , (11) 

p,a 

The following analysis is different for small and large 
values of Ec/{AE). At Ec < AE, the energy En == 
fiv{gn — Pp.) is independent of the level population and 
thus of the spin state. At relatively high temperatures, 
F < T < AE, using Eq. ([TT|), one obtains 



/ 



2 2^ 



tanh ■ 



El 
2T 



tanh ■ 



E^ 
2T 



, /o . iirp. (12) 



En = En ± eV/2, n runs over even and odd states. The 
IV-cuives are shown in Fig. [2l The current exhibits steps 
at eV = 2{AEn ± SE) where SE = {hv/TT)[pf, - irM/d] 
is the deviation from resonance between the zero- voltage 
Fermi energy in the leads and one of the levels, M being 
an integer. For T <^ 6E and small voltages, eV < SE, 
the current is zero. For the resonance, SE = 0. the 
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FIG. 2: (Color online) (a) Current steps for zero Coulomb 
energy and T = 0. E„ — AEn is the resonance energy, SE = 
hva' /d is the deviation from the resonance, (b) For exact 
resonance, SE = 0, the steps are doubled. 



current is given by the term with n = M in the sum, 
/ = /otanh(ey/4r). At eF < T the conductance is 
Ohmic, G = cIq/^T, which agrees with Rcfs. [i,[l3. 

For large Ec, the states within the interval 5p ~ \T\'^ /d 
near the lowest resonance level, i.e., those which con- 
tribute to the current, have energies ~ Ec- Therefore 
the distribution function sign((7„ —p^j.) can be used up to 
temperatures T ^ Ec, regardless of the relation between 
T and T or AE. Using Eq. ^ we find 



I={Io/i)J2mPn-p^)-HPn 



-P',)] 



(13) 



The level positions are determined by the condition of 
minimal total energy which fixes the value of iV — iVo. 
Thus the levels, as functions of the bias voltage, will cross 
the Fermi energy pairwise, one level from above while the 
other from below, keeping the number of electrons N—Nq 
unchanged and lifting the Coulomb blockade in the bias 
voltage. If the Fermi level lies between the resonances, 
the first step of hight Iq in the current appears for V = 
AE/e. The next steps appear when the bias voltage is 
increased by SV ~ 2AEn/e. The heights of the current 
steps arc the same as for low Ec'- the extra factor ^ 
in Eq. (|13p is compensated by the pairwise level crossing. 
However, the fine structure shown in Fig.l^Ja) disappears. 
If the system is close to the degeneracy point N — 
[No] = ±1/2, when one of the values p„ is close to Pf_,, 
even a small bias voltage is sufficient to produce current. 
We find from Eq. (fT3|) for low voltage, eV <C F, 



GoV 
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2 |T|4 + 4(p„-p,0'rf'' 



Go 
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(14) 



In the resonance, p„ = p^, the conductance is half of 
the conductance quantum, Go/2. The conductance does 
not depend on temperature at T < Ec- This fact differs 
from the tunnel-approximation result of Ref. [6| for the 
temperature domain T <S^ T <S^ Ec- The difference is 
due to the step in the energy spectrum caused by rela- 
tively strong coupling to the leads. Probing the onset of 
Ohmic conductance at the degeneracy point allows one 
to monitor the effective number of electrons in the QC. 



Conclusion. We develop a mean field type description 
of the Coulomb effects in the "weakly open" ID systems 
and analyze both the excitation energy spectrum and the 
electric current. The IV curves show a step- like depen- 
dence on the bias voltage, the exact shape being deter- 
mined by the ratio between the charging energy and the 
level spacing. At large charging energy and low temper- 
atures T <^ Ec, the low- voltage Ohmic conductance, 
Eq. p4|) . is temperature-independent irrespectively to 
the magnitude of the ratios T/ AE and T/T. 
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